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INTIODUCTION 


There  has  been  considerable  inte’^est  recently  in  obtaining  accurate 
numerical  values  for  the  stress  Intensity  factors  in  the  vicinity  of 
crack  tips  In  tensile  sheet.  Previous  to  several  recent  solutions, 
investigators  In  the  field  of  fracture  nechsmlcs  relied  considerably  on 
Irwin's^  approxlmati jn  derived  from  Westergaard's*  solution  for  a  series 
of  equally  spaced  collnear  creeks  in  an  infinite  sheet. 

Recently,  one  of  the  authors  has  succeeded  ir  obtaining  accurate 
results  for  the  case  of  symsctric  edge  cracks  in  rectanpilar  tensile  sheet 
by  using  a  complex  variable  approach  with  a  complex  mapping  function  to 
describe  the  gecnetry.  The  flexibility  of  the  appi^oach  was  Illustrated 
by  the  relative  simplicity  in  which  the  length/width  ratio  parameter  for 
finite  rectem^ilaLT  sheet  was  included  in  the  analysis-  The  success  of  the 
analysis  relied  heavily  on  eui  effective  plan  invented  for  the  truncation 
of  the  series  representation  of  the  exact  mapping  function. 

In  this  report,  the  technique  is  applied  to  the  problem  of  a  single 
edge  crack  in  rectangulsx  tensile  sheet.  The  analysis  is  complicated 
somewhat  by  a  certain  loss  of  mathematical  symmetry  as  'ompared  with  the 
previous  work.  For  deep  cracks  the  numerical  results  are  considerably 
different  from  the  previous  cases  thereby  indicating  the  sensitivity  of 
the  stress  intensity  factor  to  bending. 

INITIAL  FONMULATION 

The  rectangular  sheet  under  tension  weakened  by  a  single  edge  crack 
(Flgire  l)  will  be  considered  as  lying  in  the  complex  Z-plane,  Z  =  x  ly, 
with  the  center  of  the  sheet  described  by  Z  =  Z^  where  is  real . 

The  complex  variable  methods  of  Muskhellshvlll®  depend  on  the  repre¬ 
sentation  of  the  well-known  Airy  stress  function  U(x,y)  in  terms  of  two 
analytic  iunctions  of  the  complex  variable  Z  namely,  tp(z)  and  Y(z),  where 

U(x,y)  -  Re[Z<P(z)+f^(z)dZ].  (-) 
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To  facilitate  the  consideration  of  boundary  conditions,  an  eiuxiliary 
complex  plane,  the  C -plane,  is  Introduced  and  a  functional  relationship 

Z  -  u)(C)  (2) 

1 6 

is  '-und  such  that  the  unit  circle,  C  =  o  =  e  ,  and  its  interior  in  the 
C-p-cane  map  into  the  boundary  and  interior,  respectively,  of  the  region 
in  Fi»ire  1.  The  mapping  function  u)(C)  is  analytic  interior  to  the  unit 
circle  but  contains  singilarltles  necessary  to  the  description  of  corner 
points  on  the  unit  circle  itself. 

The  stress  functions  9(z)  euid  i'(z)  can  be  considered  as  function*, 
of  C*  New  notation  can  be  minimized  by  designating  5p(z)  s  <p[«-'(C)]  as  (p(C)) 
etc.,  which  leads  to  such  definitions  as  =  'p'(C)/u)'(C),  etc.  (Primes 

are  used  to  denote  differentiation.)  Thus,  the  stresses  and  displacements 
in  rectangular  coordinates  can  be  written  as 

Cy  +  0^  =  4  Re[!p'(C)/u)'(C)]  ,  (3) 

V  ■  °x  *  [v'c :)/*«'( c)]'  t  Y'(c)]V(c)  (i*) 

24  (utiv)  -  MC)  -  u<c)  ?Tcy  (^) 

where  4,  T]  are  constants  depending  on  the  material  and  bars  denote  complex 
conjugates • 

It  is  the  primary  concern  of  this  paper  to  detennlne  the  stress 
intensity  factor  at  the  crack  tip.  The  stress  Intensity  factor  K  is 
defined  by 

‘’max  "  “v  ~  t/2)(sln  3t/2)]  (6) 

'  /5? 

where  (r,t)  are  polar  coordinates  defined  with  the  origin  at  the  crack  tip 
It  can  be  shown®  that  in  terms  of  the  present  formulation. 
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K  =  2  <p'(l)/  • 


(7) 


ThUB,  the  major  effort  will  involve  determining  the  stress  function  cp(C) 
and  finding  the  numerical  values  of  !p‘(l). 

THE  MAPHlNfi  FUNCTION 

By  application  of  the  Schwartz -Christoffel  transformation,  it  was 
found  that  the  required  mapping  function  corresponds  to  an  appropriate 
branch  of 


z  -  <:)  = 


_ (C-l)  dj _ 

y(C*-2Cco6  B+1)(C®-2Cco6  a+l)(C’+l) 


(8) 


In  Equate  (.  8,  0  <  p  <  a  <  11/2,  where  a  and  P  can  be  varied  to  obtain 

desired  xatios  and  L/x^.  The  choice  of  branch  was  made  by  defining 

o/l)  =  Z  -  +  L,  thus 

'  '  o  o 

=  ^o  ■  *0  •  ^^o'  =  ^o  *  *0 

a)(e^"/2)  =  -  ly^,  etc. 

By  considering  the  mapping  function  on  the  unit  circle,  the  length 
parameters  can  be  expressed  in  terms  of  standard  elliptic  Integrals  of  the 
first  kind,  F((P,K)  [Seference  7],  as  follows: 


^0^2  ^  (l+cos  a)]  F(ii/2,  kj). 

(9) 

y^  [cos  p(l+cos  r(ll/2,  k,  ), 

(10) 

-I/2 

L  =  [cos  p  (l+cos  a)l  F(<Pj  ,  1^1)  » 

(11) 

where 
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qii  -  3ln’^  [(l+coB  a)(l-co8  p)/(l-coe  a)(l-fcoB 
kj  «  [cob  a  (1+COB  P)/cOB  P  (1+COB 
k,  =  [(cob  P-cob  a)/coB  p(lit:0B 
It  iB  convenient  to  define  <P,  by  the  relation 
9,  «  Bln‘^1  . 

Then  It  followB  ^JBnedlately  that 

kf  •  COB  . 

The  length/width  ratio  of  the  plate  now  can  be  written  aa 

ypAp  -  2  ZL»ln  ^  ,  t/2)  .  (12) 

F(eOB  H/2)  ^ 

It  iB  clear  fron  (12)  that  a  fixed  y^/x^  rat^o  detemlneB  <P,  and  thUB  a 
fixed  relation  between  a  and  p.  miB  ccnputatlon  can  be  carried  out  by 
Bl*ply  interpolating  fran  exletlng  tablee.’  Variation  of  p,  say,  then 
pendta  a  atudy  of  varloue  ration. 

The  napping  function  haa  bIx  branch  polntB  falling  on  the  unit  circle 
(Flgire  2).  In  addition  to  tbeae  comer -deecrlblrg  slngilarltleB,  the 
crack  tip  1b  deacrlbed  by  the  root  of  ti»'(o)  =  o  which  occutb  at  ct  =  1. 

The  mapping  function  can  be  expressed  In  series  form  as 

^  1  «°  (13) 
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where  the  mapping  coefficients  are  real.  Although  the  A^'s  can  be  related 
to  integrals  involving  the  elliptic  functions,  it  is  simpler  cooiutationallor 
to  calculate  them  frcn  easily  derived  recursive  f omulae . 

It  was  found  that  retention  of  the  exact  mapping  function  and  power 
series  representation  of  the  stress  functions  leads  to  the  necessity  of 
solving  an  infinite  system  of  linear  equations.  The  truncation  procedure 
developed®  was  used  in  order  to  achieve  an  effective  rate  of  convergence. 

In  general,  the  procedure  involves  the  use  of  selected  truncations  of 
the  series  (I3)  preserving  the  properties  of  the  first  and  second  deriv¬ 
atives,  ii)'(C)  and  ii)'(C),  at  the  crack  tip,  namely  at  C  =  1- 
In  peirtlcular,  we  define 


where 


•n  "  '^n’  "  ' 


m; 


(11*) 


(15) 


and 


where 


*m+l  "  *m+2  °  ® 


(m+2)S  =  Q  - 

(m+l)R  -  -Q  +  -  (m+1) 


m  ^  \ 

T  =  I  n  e  ,  =  E  n(n-l)  t 


(16) 


m  n=l  ^  ® 

Q  -  U.'(l)  =  (1/2)  ( 


/2(1-co8  P)(1-cob  a) 


). 


With  this  definition  of  u)j(C),  satisfaction  of  '*ij'(l)  =  0  and  u)  ^(l)  =Q 
is  assured.  The  choice  of  it  is  made  from  an  examination  of  the  partial  sums 
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derived  from  the  exact  mapping  function  and  selecting  those  values  of  m 
for  which  u)'^)  0  and  is  (l)  rs  Q  simulLaneously .  Thus,  the  necessary 

corrections  of  the  geometry  In  the  vicinity  of  the  crack  tip  can  be  made 
with  little  disturbance  of  the  ove*-  -all  configuration . 

DETERMINATION  OF  THE  STRESS  FUNCTIONS 

The  loading  conditions  In  Figure  1  can  be  expressed  conveniently  In 
terms  of  the  force  resultant .  Considering  an  arc  of  the  material  with 
element  dS,  we  denote  the  horizontal  and  vertical  forces  as  XdS  and  YdS, 
respectively.  If  the  arc  Is  taken  as  the  boundeuy  of  Figure  1,  then  S  can 
be  considered  as  a  function  of  o.  The  boundary  condition  can  be  written  as 

cp(o)  +  u)(o)  cp'(o)A'(o)  +  T(a)  =  J®(x+lY)dS  =  g(o).  (17) 

The  solution  then  requires  the  determination  of  the  functions 
<p(C)  and  tCC)  which  are  analytic  for  |  Ci<l  and  satisfy  the  loading 
conditions  (Equation  I7). 

The  problem  will  be  considered  in  terms  of  polynomial  approximations 
of  the  exact  geometry  as  Indicated  by  the  truncation  procedure  of  the 
previous  paragraph . 

It  was  shown  In  Reference  3  that  the  stress  functions  cp(C)i  Y(C)  have 
the  following  general  structures: 

cp(t)  =  (tA)  + 'Pi(C)  (18) 

T(C)  =  (T/2)  a)(C)  +  Y,(C)  (19) 

where 

P  a>  „  m+2  -1 

91(0=1  E  c  C  +  E  a  C" 

[n^  n  n  J 

'^l(C)  =  t[a/(c-i)  + 

L  n=0  “  J 

and 

— • 
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with  the  above  forme  of  9(C),  buindary  condition  (ij)  may  be 

vrritten  as  follows; 


9j(a)  +  u)(o)  9j'(o)  /  is'(o)  +  Y^Co)  =  g(a)  -  T/2[u)(a)+u)(o)]  =  G(o), 

(20) 


If  G(cj)  is  expanded  in  the  Fourier  series, 

G(a)  =  ?  C„  0^  (21) 

K=-od  ^ 

then, 

Cj^=  ^[aj(e^^)+uo(e^^)]  cos  nej(co6  KG-l  sin  Ke)d^ 

^  ^-K  for  K  =  1,  2, . .  . 


Insertion  of  the  appropriate  series  expansions  into  Equation  20  and  equating 
coefficients  of  equal  powers  of  a  (from  the  positive  powers),  one  finds  the 
system  for  determination  of  to  be 


o+i-P  ,  , 

b  n(€Ci!.,+C  •j_t)=0, 

n  =  l  '  n  n+p-1  n  n+p-1  n  n+p-1' 


P  =  1,  2, 


. ,  m+2 . 


(22) 


The  calculation  of  the  stress  intensity  factor  K  requires  the  determination 
of  9'(l)  as  Indicated  in  Equation  7 

From  Equation  l8  it  can  be  seen  tiiaL 


'P'(l) 


mt2 
+  E 
n  =  l 


=T  [Ej  t  5^]  . 


(23) 


Ej  can  be  determined  directly  from  the  Fourier  coefficients  of  the 
applied  load.  5^  can  be  calculated  from  the  solution  of  the  system  of 
Equations 
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The  analyels  was  carried  out  numerically  for  the  two  cases  y  /x  =1.0 
^0^*0  "3-06.  Ir.  general,  three  values  of  m  were  selected  to  study  the 
convergence  of  Ej  and  .  The  range  of  m  considered  fell  In  the  Intervals 
0  <  mi  <  50,  50  <  m,  <  100,  100  <  m,  <  150.  A  study  of  the  convergence  of 
Ej  and  led  to  the  assertion  that  the  calculated  values  of  ip*(l)/T  su'e 
numerically  correct  to  within  2  percent . 

The  key  data  eu^  presented  In  Table  I  below: 


TABLE  I 

BVALUATIOH  OF  K/T 


a 

P 

X 

0 

0 

X 

0 

>» 

L/2x 

0 

tp'Cl)! 

k/t 

K/K, 

O.3706 

0.1600 

1.150 

0.997 

0.11*1* 

1.19 

0.69 

1.19 

0.41*  1*5 

0.3000 

1.178 

0.997 

0.254 

1.33 

1.15 

1.45 

0.7573 

0.7000 

1.381 

0.997 

0.487 

1.70 

2.88 

2.22 

1.0279 

1.0000 

1.75'* 

0.997 

0.621 

2.49 

5.75 

3.17 

1.21^1 

1.2000 

2.271 

0.997 

0.711 

3.74 

10,11 

4,16 

1.5095 

0.5000 

0.81*1 

3.056 

0.170 

0.331 

0.64 

1.18 

1-5255 

1.0000 

1.080 

3.056 

0.378 

0.604 

1.65 

1-71 

1.5368 

1.2000 

1.326 

3.056 

0.490 

0.95'* 

2.84 

2.22 

1-5525 

1.1*000 

1.951 

3.056 

0.647 

2.154 

6.88 

3.44 

The  final  column  Is  presented  for  the  seike  of  comparison  with 

Irwin's  approximation  for  symmetric  edge  cracks.  In  the  present  notation. 


f  „  1 

V'  ■  I  IT  ■  (2'*) 

The  ratio  K/k^  as  a  function  of  L/2x^  Is  shown  in  Figure  3.  Its 
departure  from  unity  is  dcralnantly  due  to  the  bending  effect  of  the  applied 
load. 
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DISCUSSION 


For  very  short  crack  lengths  one  would  expect  the  solution  to  behave 

as  if  the  plate  were  semi -inf ini te  in  extent.  This  is  consistent  with  the 

numerical  result  K/K.,  1.13  as  L  -•  0  in  Figure  3- 

A 

Comparison  of  the  results  for  Fq/Xq  *  1.00  and  Xq/Xq  '  3-06  indicates 
agreement  for  very  short  and  fairly  long  crack  depths.  The  variation  in 
the  interval  0.2  <  L/2Xq  <  0.5  it  due  to  the  difference  in  length/width 
ratios  .  It  is  reasonable  to  expect  that  the  results  would  not  change  for 
larger  F.^/Xq  ratios  . 
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